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Abstract. Orbital period changes are an important diagnostic for understanding low mass X-ray binary (LMXB) accretion- 
induced angular momentum exchange and overall system evolution. We present our most recent results for the eclipse 
timing of the LMXB EXO 0748—676. Since its discovery in 1985 it has apparently undergone three distinct orbital period 
"epochs", each characterized by a different orbital period than the previous epoch. We outline the orbital period behavior for 
EXO 0748— 676 over the past 18 years and discuss the implications of this behavior in light of current theoretical ideas for 
LMXB evolution. 



INTRODUCTION 

There are eight currently known low mass X-ray binary 
(LMXB) systems that undergo full or partial eclipses. 
Such systems are important for the study of evolution in 
LMXBs because the eclipse edges provide timing mark- 
ers that make possible the systematic observation of or- 
bital period changes. The best studied orbit in this group 
is that of EXO 0748— 676 which has been X-ray active 
since 1985 [1]. This long timeline has allowed an un- 
precedented look at its orbital dynamics [see 2, and ref- 
erences therein]. The emerging picture for the orbital 
period behavior, however, is anything but the expected 
smooth variations in P or t> and P or b based on theoreti- 
cal calculations done to date [e.g., see 3, and references 
therein]. Rather, the observed period changes are discon- 
tinuous across multiple distinct epochs, and large appar- 
ent changes in P ot ^ of the order seconds can be observed 
on timescales as short as one orbit. The magnitude of the 
observed changes in orbital period are much larger than 
expected from LMXB evolutionary theory. This is likely 
an indication that the observed variations in P or b are short 
timescale effects of angular momentum redistribution in 
the system and are masking the underlying long term or- 
bit evolution. The Rossi X-Ray Timing Explorer (RXTE) 
satellite has allowed us to monitor systematically the or- 
bit of EXO 0748—676 in an effort to delineate and under- 
stand the process of angular momentum exchange and its 
effects on the EXO 0748— 676 orbit. In this conference 
paper we report our recent progress. 



RXTE OBSERVATIONS OF 
EXO 0748-676 

We now have 85 additional RXTE mid-eclipse time mea- 
surements over and above those analyzed in Wolff et al. 
[2] for EXO 0748 -676. Wolff et al. reported that the 
orbit period increased around 1990 in a somewhat dis- 
continuous manner by ~7.9 ms since its discovery in 
February 1985. The cumulative data set of mid-eclipse 
timings considered in that paper ended in December 
2000. Our present data set extends to observations done 
through October 5, 2003. O-C (Observed - Calcu- 
lated) diagrams for all published mid-eclipse timings of 
EXO 0748— 676 plus the additional timings we present 
here (a total of 262 mid-eclipse time measurements) are 
shown in Figures 1 and 2. 

The additional eclipse measurements display a re- 
markable behavior. When the eclipses observed after De- 
cember 2000 are included in our analysis the orbit pe- 
riod appears to return to a value close to the orbit period 
observed during the EXOSAT/GINGA epoch. We fit a 
three-constant-period model (P O rb0, ^orbi, and P rb2) to 
these data in which we constrain the phase to be constant 
across the instantaneous period changes but let the cy- 
cles («£o and n^i) °f the period change be free parameters 
(see Table 1). The best-fit model gives a positive period 
change of about P or t>i — P O rb0 = 8.19 ± 0.06 ms similar 
to that reported in [2] during the GINGA to early RXTE 
era. However, the additional data does not continue the 
trend of random wandering about the new period. In- 
stead, the apparent EXO 0748—676 orbital period under- 
goes another change nearMJD 51710 (near cycle 35136) 
of ^01*2 — ^orbi = —6.75 ±0.02 ms. This brings the 
net change in orbit period in the EXOSAT— RXTE eras 
(February 1985 through October 2003) to only P G rb2 - 
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FIGURE 1. The O — C residuals for all published mid-eclipse timings of EXO 0748— 676 plus those additional timings we 
present here for a total of 262 mid-eclipse timings. All times are in Modified Julian Days corrected to the solar system barycenter. 
The measurement errors are as shown and except for a few cases are less than ~ 1 — 2 seconds. The solid line represents the Double- 
Broken Constant Period model for these data discussed in the text. The large deflection of the ASCA data away from either solution 
may be attributable to large excursions due to accumulated jitter as discussed in [2]. The systematic wandering of the orbit period 
caused by cumulative jitter prior to MJD 51710 is apparent. After this date, however, the orbit period takes on a new value that is 
only 1 .4 ms larger than the apparent orbit period prior to MJD 48026. 
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FIGURE 2. The O - C residuals for all USA and RXTE mid-eclipse timings of EXO 0748-676 plotted against the middle 
period of the Double-Broken Constant Period solution in Table 1. All times are in Modified Julian Days corrected to the solar 
system barycenter. The solid line represents a constant period solution for the O — C residuals and is clearly unacceptable as a 
model for these data. The systematic wandering of the orbit period caused by cumulative jitter prior to MJD 51710 is apparent. 
After this date, however, the orbit period takes on a new value of P orb = 0.15933774511(16) x 10~ 10 (days) and the wandering 
about the mean period is noticeably reduced. 
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P or t,0 = 1 .45 ±0.06 ms, significantly less than reported in 
[2]. Note that the orbit period in the present epoch (P rb2) 
is distinct from the period in the EXOS AT/GINGA epoch 
(^orbo) according to our analysis. 



CUMULATIVE JITTER IN X-RAY 
MID-ECLIPSE TIMINGS 

The O C residuals in Figure 2 appear correlated, that is, 
they are often systematically grouped above and below 
the zero residual line. In some cases the positive or nega- 
tive variations in the mid-eclipse residuals can persist for 
many months. We have referred to this variation in the 
apparent orbital period as "jitter", and apparent changes 
in orbital period can be caused by accumulated intrinsic 
jitter in a constant orbital period. "Measurement error" 
is simply the random error associated with our measure- 
ment of the individual mid-eclipse times. We take "in- 
trinsic jitter" to mean the systematic wandering of the 
mid-eclipse timings around a smooth underlying model 
that can not be accounted for by random uncorrected 
measurement errors. In particular, the systematic wan- 
dering of the mid-eclipse residuals apparent in Figure 2 
can be accounted for by the accumulated intrinsic period 
jitter in individual orbit cycles. In order to investigate 
the observational measurement error inherent in our tim- 
ing analysis and the intrinsic period jitter that is caused 
by other mechanisms (possibly mechanisms inside the 
LMXB system) we apply Maximum Likelihood Method 
[MLM, see 4] to estimate the parameters of a model 
for the orbital evolution that simultaneously includes an 
orbit period, an orbit period derivative, non-zero ran- 
dom intrinsic scatter (variance a|), and non-zero random 
measurement error (variance of). We give results for 
EXO 0748-676 in two cases: (a) RXTE and USA data 
prior to MJD 51710, and, (b) the RXTE data after MJD 
51710. A MLM model for the EXO 0748 -676 RXTE 
and USA mid-eclipse timing data prior to MJD 51710 
that includes a non-zero period derivative, intrinsic jitter 
and measurement error is not statistically preferred over 
a model with only intrinsic jitter and measurement er- 
ror (a model with no intrinsic period jitter is strongly re- 
jected according to the MLM analysis). The MLM anal- 
ysis suggests that the intrinsic jitter in the eclipse timings 
is characterized by a e ~ 0.13 s with measurement error 
a e ~ 1.3 s, similar to the results in [2]. However, for the 
RXTE data after MJD 51710 the MLM analysis suggests 
that while the measurement error in the eclipse timings 
remains roughly constant at a e ~ 1 .2 s, the intrinsic pe- 
riod jitter is greatly reduced to a £ ~ 0.028 s. Also, as 
was the case for the pre-MJD 51710 data, a model for the 
post-MJD 51710 timing data that includes a non-zero pe- 
riod derivative, intrinsic jitter and measurement error is 



not statistically preferred over a model with only intrin- 
sic period jitter and measurement error. Whatever pro- 
cess resulted in the abrupt change in the orbital period 
after MJD 51710 also brought about a reduction in the 
magnitude of the intrinsic orbital period jitter according 
to the MLM analysis. 



DISCUSSION 

The X-ray eclipse timings for EXO 0748— 676 appear 
to show three distinct orbital periods in three succes- 
sive epochs along with significant intrinsic period jit- 
ter in these epochs. Furthermore, inspection of Fig- 
ure 2 shows that during the early part of the RXTE 
epoch (until approximately March 2000) the period jit- 
ter is especially prominent. After this date the orbit pe- 
riod changes again by ~-6.7 ms nearly returning to the 
value characteristic of the EXOSAT-GINGA epoch. 
Such abruptly changing orbit periods are observed in 
a number of Algol binaries [see the discussion in [5] 
and references found therein] although with larger ampli- 
tude in AP or b/.Porb. F° r Algol systems AP or b/P rb ~ 10~ 5 
whereas in EXO0748-676 we find AP orb /P orb - 10~ 6 . 
[6] attributed the abrupt period changes in Algol systems 
to magnetic activity in a convective secondary star in- 
ducing changes in its quadrupole moment and altering 
the orbital angular momentum. If this explanation is ap- 
plicable to the period changes in EXO 0748— 676 [see 
also 7] then around the March 2000 magnetic activity 
in the secondary was altered in some manner changing 
the secondary's influence on the system orbital angular 
momentum distribution. Before this date the orbit period 
jittered around a mean value but the jitter is reduced af- 
ter this date. This may imply that magnetic activity in 
the secondary also was reduced, perhaps as a result of 
magnetic cycling similar to the 1 1 -year cycle in the Sun. 
Small changes in radius of the secondary are predicted by 
the magnetic activity explanation of Algol orbit period 
changes put forward by Applegate and Patterson [8] due 
to the changing magnetic pressure support for the outer 
layers of the secondary star. The data in Figure 3 suggest 
that the radius of the occulting star decreased across the 
MJD 51710 boundary. Totality duration is ^7 seconds 
shorter after MJD 51710 than before MJD 51710. How- 
ever, a spherically symmetric decrease of 1-2% in the 
radius of the secondary star is unlikely and other non- 
spherical changes in the radius of the secondary must be 
considered as well. 

The measurements shown in Figures 1, 2, and 3 give 
a magnitude and direction for the abrupt orbital period 
changes and the magnitude and direction of the abrupt 
change in the duration of X-ray totality. Any attempt 
to understand these data will have to pose some sort 
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FIGURE 3. The durations of all RXTE observed eclipses plotted as a function of observed O C residual for EXO 0748-676. 
The points at increasing negative O - C values are the same eclipses observed after the dramatic change in period near MJD 51710 
in Figure 2. 



of almost instantaneous change in the binary system 
parameters with an accompanying redistribution of the 
orbital and spin angular momenta rather than a slow 
change in the system because of the effects of mass 
exchange. 
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TABLE 1. Orbital Ephemerides 
EXO 0748-676 



of 



Parameter 



Value 



All Eclipse Timing Data 



T Q (MJD/TDB) 


= 46111.0752010(42) 


Jorbo (day) 


= 0.15933772838(66) 


«m (cycle) 


12019.5(63.8) 


^orbi (day) 


= 0.15933782322(10) 


n bl (cycle) 


= 35136.4(11.4) 


Porbl (day) 


= 0.15933774511(16) 


X 2 (dof) 


46.1(256) 
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